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1. INTRODUCTION 

In this work minimizing true power loss is the main objective of the problem. A variety of methods 
[1-6] have been applied to solve the problem. Subsequently various evolutionary methods [7-16] applied to 
solve the problem, in that many algorithms stuck in local optimal solution In this work amplified brain storm 
optimization (ABS) algorithm and quantum based brain storm (QBS) optimization algorithm is used for 
solving optimal reactive power problem. Brain storm optimization (BSO) algorithm gets trapped into local 
optima when applied to different optimization problems. In the mathematical field of graph theory, a 
Hamiltonian path is a path in an undirected or directed graph that visits each vertex exactly once. In the 
proposed algorithm Hamiltonian cycle will improve the explore abilities and also stay away from local 
optimal solution. In QBS algorithm completely, the mechanism of quantum behavior, which causes uncertain 
of every idea lead to a superior capability to bounce out of the local optimal solution. Proposed ABS 
algorithm and QBS optimization algorithm has been tested in standard IEEE 57 bus test system. 


2. PROBLEM FORMULATION 
2.1. Real power loss 


F = P, = Yxenpr 8k (V? + VŽ — 2V; Vjcos8;;) (1) 
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2.2. Amplification of voltage profile 
F = P, + œ X Voltage Deviation (2) 
Voltage deviation given by: 
Voltage Deviation = bye PA; -1| (3) 
2.3. Constraint (equality) 
Pg = Pp +P, (4) 


2.4. Constraints (inequality) 


Pence = Pgslack Pease (5) 
On” < Qai QP ie Ng (6) 
ymin < y; < VMA j €N (7) 
Te < T; < TP? i € Ny (8) 
QE" < Qe < QE*,1 E Ne (9) 


3. AMPLIFIED BRAIN STORM OPTIMIZATION ALGORITHM 

BSO [17] gets trapped into local optima when applied to different optimization problems. In the 
projected amplified brain storm optimization algorithm Hamiltonian cycle has been applied to improve the 
search abilities and also to avoid of trap in local optimal solution. A node is arbitrarily chosen from the graph 
as the preliminary point to form a Hamiltonian cycle. At generation t and t+1, L; and Lı are the length of 
Hamiltonian cycle correspondingly. Their ratio r at generation t(r,) can be described as: 


Lt (10) 


C= 
t Le 


Hamilton cycle algorithm as follows: 


Commence 

Step 1: node vl chosen as initial point,. 

Step 2: V4,V2,..,V; is chosen and V4, is picked with least weight linkingv,, then the 14,0 ,..,Vj41 
is obtained. 

Step 3: when it+l<n, subsequently itl is used to substitute i, and revisit to Step 2; 
condition not occurred , then revisit to the final Hamiltonian cycle C=1,,03,,..,UV,¥,then go 
back to Step 4. 

Step 4: vv, EEG) ; Viz Vj41 €E ECG), w(viv;) + WhViVj1) < W(ViVi+1) + w(v;%;41) 

Then C, = (C — {vjvi41, VjVj+1}) U (vir), Vierdjaa} 

End if 

End for 

Step 5: C is substituted by Cl, and revisit Step 4. 

Step 6: compute the extent of the Hamiltonian cycle C. 

End for i 


In the proposed amplified brain storm optimization (ABSO) algorithm Hamiltonian cycle will 
improve the explore abilities and also stay away from local optimal solution. 


Commence 

Step 1: “n” potential solutions are arbitrarily engendered 

Step 2: “n” individuals are clustered into “m” clusters 

Step 3: “n” individuals will be appraised 

Step 4: In every cluster rank the individuals then the most excellent individual’s are 
recorded as cluster center Step 5: Between 0 and 1 arbitrarily a value will be engendered; 
If the value is smaller than a probability; then 

i. a cluster center has been Arbitrarily chosen; ii. To swap the certain cluster center 
arbitrarily engender an individual 
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Step 6: new-fangled individuals are engendered 

Calculate the Hamiltonian cycle C and its extent Li by Hamilton algorithm 

{ 

Commence 

Step 1: node vl chosen as initial point, 

Step 2: Vi V2.. Vi is chosen and Vi} is picked with least weight linkingv,, then the 
V1,V2,..,Vig4, is obtained. 

Step 3: when itl<n, subsequently i+1 is used to substitute i, and revisit to Step 2 
Step 4: for all i and j in cycle C, if 1< i+l<j<n, then 

i +j; viv; € E(G) 

Visa Vjar E E(G),w(viv;) + w(VirYj) < w(ViYVi) + w(v;V;41) 

Then C, = (C — {vjvi41, VjVj+1}) U (vir), Vierdjaa} 

End if 

End for 

Step 5: C is substituted by Cl, and revisit Step 4. 

Step 6: compute the extent of the Hamiltonian cycle C. 

End for “i” 


Å 
When t>1 then calculate value of the ņ by n= 
¢ 


End if 
Execute decision optimization procedure 
{ 


Commence 


Tminimum < Tt < Tmaximum OTF Tt = maximum 
Arbitrarily engender n, individuals; 
End if 
End 
} 
Calculate the population according to the recently modernized positions; 
t = t+1. 
Step 7: when “n” new-fangled individuals are engendered, then go to Step 8; or else go to 
Step 6. 
Step 8:end conditions met ; or else go to Step 2. 
End 


4. QUANTUM BASED BRAIN STORM OPTIMIZATION ALGORITHM 

In BSO algorithm population is indicated as swarm moreover every individual is described as an 
idea. Originally, every idea is arbitrarily initialized inside the exploration space. Subsequently most excellent 
one in every cluster is selected as the cluster centre. Sporadically, an arbitrarily chosen centre is swapped by 
a recently engendered idea, by that the swarm has been kept away from the local optimum. 


ij 
new 


X x, +E N(u,0) (11) 


ij 
new 


xj xij 
x. W1Xo1a1 + W2Xo1a2 (12) 


č is a factor used in the evolution process and can be articulated as, 


(13) 


(mAN >j random 


&(t) = logsig 
Quantum state of an idea is illustrated by a wave function ¥ (x, t) as an alternative of the position 
modernized only in Brain storm optimization algorithm. By using Schrödinger equation probability density 
function of the position is identified such that each idea is located. Monte Carlo simulation method is used, to 
measure the position for each idea from the quantum state to the traditional one. 


y 4 + (14/2) n/n) (random < 0.5) 


new lay — (1/2) m/w) (random = 0.5) (14) 
qij = random*x) peste + A - random)*x!),.., (15) 
lj = 2b|mbest/ a zdal (16) 
b= 1-05 (17) 
mbest) = k axt o/k (18) 
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random * xf ost + (1 —random) * A hog + (bIZE, xË ce /k - xl) * In(1/u) + € N(u,o) (random < 0.5) (19) 


random * xf ost + (1 — random) « gos + (blz, x cel k - xl) x In(1/u) + ë N(u, o) (random = 0.5) 


U 
Xnew = 


Step a: Initialize the parameters. 

Step b: Arbitrarily produce “n” ideas Step c: By k-means algorithm cluster ideas. 
Step d: With a predetermined probability modernize the centre of a capriciously chosen 
cluster. 
Step e: Individual generation created. 

Step f: Quantum mechanism is exploited based on the chosen idea 

Step g: Crossover operator is implemented 

Step h: evaluate the new-fangled idea with the older one, 

Step i: If “n” ideas have been engender, then go to Step 9. Or else go to Step 5. 

Step j: Stop whether the present number of iterations Nc attain the Nmax. or else, go to 


wa 


n 


5. SIMULATION STUDY 
Proposed ABS optimization algorithm and QBS optimization algorithm has been tested, in IEEE 57 
Bus system [18]. Table 1 shows the comparison results. 


Table 1. Simulation results of IEEE-57 system 


Control variables Base case MPSO [19] __ PSO[19] _CGA[19]_ _AGA[19]__ ABS QBS 
Val 1.040 1.093 1.083 0.968 1.027 1.019 1.020 
VG2 1.010 1.086 1.071 1.049 1.011 1.025 1.022 
VG3 0.985 1.056 1.055 1.056 1.033 1.027 1.019 
VG6 0.980 1.038 1.036 0.987 1.001 1.021 1.012 
VG8 1.005 1.066 1.059 1.022 1.051 1.027 1.037 
VG9 0.980 1.054 1.048 0.991 1.051 1.035 1.028 
VG 12 1.015 1.054 1.046 1.004 1.057 1.049 1.046 
Tap 19 0.970 0.975 0.987 0.920 1.030 0.908 0.900 
Tap 20 0.978 0.982 0.983 0.920 1.020 0.906 0.911 
Tap 31 1.043 0.975 0.981 0.970 1.060 0.909 0.916 
Tap 35 1.000 1.025 1.003 NR* NR* 1.013 1.014 
Tap 36 1.000 1.002 0.985 NR* NR* 1.015 1.012 
Tap 37 1.043 1.007 1.009 0.900 0.990 1.006 1.017 
Tap 41 0.967 0.994 1.007 0.910 1.100 0.947 0.936 
Tap 46 0.975 1.013 1.018 1.100 0.980 1.019 1.014 
Tap 54 0.955 0.988 0.986 0.940 1.010 0.921 0.920 
Tap 58 0.955 0.979 0.992 0.950 1.080 0.937 0.932 
Tap 59 0.900 0.983 0.990 1.030 0.940 0.926 0.921 
Tap 65 0.930 1.015 0.997 1.090 0.950 1.006 1.013 
Tap 66 0.895 0.975 0.984 0.900 1.050 0.934 0.926 
Tap 71 0.958 1.020 0.990 0.900 0.950 1.006 1.052 
Tap 73 0.958 1.001 0.988 1.000 1.010 1.013 1.007 
Tap 76 0.980 0.979 0.980 0.960 0.940 0.947 0.923 
Tap 80 0.940 1.002 1.017 1.000 1.000 1.009 1.037 
QC 18 0.1 0.179 0.131 0.084 0.016 0.150 0.147 
QC 25 0.059 0.176 0.144 0.008 0.015 0.142 0.138 
QC 53 0.063 0.141 0.162 0.053 0.038 0.127 0.121 
PG (MW) 1278.6 1274.4 1274.8 1276 1275 1272.99 1272.04 
QG (Mvar) 321.08 272.27 276.58 309.1 304.4 242.51 272.12 
Reduction in PLoss (%) 0 15.4 14.1 9.2 11.6 25:32 27.88 
Total PLoss (Mw) 27.8 23.51 23.86 25.24 24.56 20.760 20.049 


NR* - Not reported. 


6. CONCLUSION 

In this paper ABS optimization algorithm and QBS optimization algorithm successfully solved the 
optimal reactive power problem. In projetced ABS algorithm to escape BSO from local optima and to 
maintain the optimization process Hamiltonian cycle has been utilized. In the mathematical field of graph 
theory, a Hamiltonian path is a path in an undirected or directed graph that visits each vertex exactly once. In 
QBS approach by using Schrédinger equation probability density function of the position is identified such 
that each idea is located. Monte Carlo simulation method is used, to measure the position for each idea from 
the quantum state to the traditional one. Proposed ABS algorithm and QBS optimization algorithm has been 
tested in standard IEEE 57 bus test system and simulation results show the projected algorithms reduced the 
real power loss efficiently. 
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